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SUMMARY 


Mathematical  modeling  is  becoming  a  vital  tool  in  Ungineerir.g  design. 
These  models  have  been  used  to  provide  decision  support  to  designers  of 
systems  characterized  by  many  design  choices.  However,  for  the  design  of 
complex  systems,  these  models  become  too  large  and  cumbersome  to  be  of 
significant  value.  As  a  result,  research  has  progressed  in  the 
decomposition  of  the  design  problem.  / 

A  decision  support  system,  serving  as  a  computer-based /aid  for 
engineers  designing  weapon  systems,  must  integrate  "procedures  dealing  with 
designing  for  producibility  and  designing  for  supportability  with  those 
dealing  with  designing  for  performance  cost  and  schedule'^  (ULCE 
Implementation  Plan,  1987).  One  of  the  characteristics  of  large  weapon 
system  development  that  causes  difficulties  in  achieving  this  integration 
is  the  hierarchical  nature  of  the  design  process.  Typically,  a  system  is 
designed  in  terms  of  subsystems  which  are  required  to  meet  certain 
specifications.  Problems  begin  to  arise  when  these  specifications  cannot 
be  met.  Consequently,  these  design  processes  may  be  enhanced  by  a 
computer-based  system  that  supports  the  hierarchical  nature  of  design. 

This  research  examines  the  application  of  mathematical  modeling  and 
knowledge-based  system  techniques  to  the  development  of  decision  support 


systems  to  support  engineering  design. 
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This  research  was  conducted  in  support  of  Air  Force  Human  Resource 
Laboratory  (AFHRL)  involvement  in  the  development  of  a  prototype  Unified 
Life  Cycle  Engineering  (ULCE)  system.  The  intent  of  this  system  is  to 
serve  as  a  computer-based  aid  for  engineers  designing  weapon  systems,  in 
particular  integrating  procedures  to  address  producibility  and 
supportability  in  weapon  system  design. 

I  wish  to  thank  AFHRL  for  their  sponsorship  and  UES  for  their  attention 
during  the  SFRP.  Several  individuals  deserve  special  identification  for 
their  support  of  my  work  at  AFHRL.  Lt.  Col.  Joe  Coleman  was  instrumental 
in  helping  me  to  identify  an  excellent  research  direction.  Capt.  Tom  King 
provided  timely  resources  and  support.  Ma j .  Joe  Nerad,  Capt.  Maureen 
Harrington  and  Lt.  Lee  Dayton  helped  me  several  times  in  overcoming 
hardware/software  difficulties. 

I  greatly  appreciate  the  hospitality  that  I  received  from  Bert  Cream 
and  Col.  Donald  Tetmeyer.  I  would  also  like  to  thank  Richard  Lamb,  Jim 
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ENGINEERING  DESIGN  WITH  DECISION  SUPPORT: 

AN  APPLICATION  OF  GOAL  DECOMPOSITION 

1.  introduction 

Engineering  design  in  its  broadest  sense  refers  to  the 
activity  of  selecting  and  placing  materials  to  form  a 
system.  Because  of  the  inherent  complexities,  design  engineers 
have  increasingly  depended  upon  mathematical  modeling  to  help 
guide  the  design  process  [12].  However,  for  the  design  of  large 
systems,  these  models  quickly  become  too  large  and  cumbersome  to 
be  of  significant  value.  As  a  result,  research  has  progressed 
into  the  decomposition  of  the  design  problem  with  the  intant  of 
breaking  the  main  design  into  a  series  of  smaller  design  problems 
that  can  be  solved  and  then  recombined  in  the  hope  of  generating 
a  good  global  design. 

Figure  1  illustrates  a  network  model  of  the  design 
decomposition  problem.  The  node  at  the  top  level  (level  0;  of 
the  network  represents  the  system  which  is  to  be  designed.  As  we 
pass  downward  througli  the  network,  the  parts  of  the  system  are 
exploded  into  greater  detail,  until  we  reach  the  bottom  of  the 
network  where  no  further  decomposition  is  possible. 

The  network  shown  in  Figure  1  can  be  thought  of  as  actually 
containing  two  different  types  of  optimization  problems:  (1)  A 
high  level  optimization,  and  (2)  Several  lower  level 
optimizations.  The  high  level  optimization  is  equivalent  to 
answering  the  question,  "What  type  of  system  do  we  want  and  what 
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do  we  want  it  to  be  able  to  do?"  An  example  of  this  type  of 
optimization  problem,  which  is  typically  characterized  by 
multiple  objectives,  can  be  found  in  [5].  In  the  case  that 
several  designs  have  already  been  established,  Pareto  optimal 
designs  can  be  determined  by  the  method  given  in  [11]  to  reduce 
to  list  of  candidate  designs  to  some  managable  number.  The  design 
decisions  that  are  made  at  this  level  are  then  passed  downward  to 
the  next  level  (level  1)  in  the  form  of  system  specifications. 

At  the  middle  and  lowest  levels  in  the  network,  the 
objective  is  not  to  optimize  the  design  of  the  subcomponents  over 
some  set  of  objectives,  but  instead  to  meet  some  stated  level  of 
attributes  that  has  been  provideu  by  a  parent  optimization.  In 
some  cases  it  may  be  suboptimal  for  a  subcomponent  to  over- 
perform  as  well  as  under-perform.  For  example,  in  designing  an 
aircraft,  we  would  like  to  achieve  a  specified  level  of  speed  (an 
attribute  measuring  performance);  however,  optimizing  the  wing 
design  for  speed  would  cleaily  be  suboptimal  if  the  landing  gear 
of  the  aircraft  is  not  designed  to  accommodate  the  greater 
speeds . 

Of  course,  other  objectives  besides  performance  may  be  of 
interest.  Reliability,  cost,  support ab i 1 i ty  and  producibi  1  i  ty 
are  some  objectives  that  frequently  must  be  addressed  in 
designing  systems.  In  each  case,  a  stated  level  of  achievement 
is  what  is  required  at  the  subsystem  level. 

The  need  for  communication  between  levels  in  the  hierarchy 
of  Figure  1  is  currently  satisfied  by  multiple  design  iterations 
in  real  world  design  processes.  In  fact,  it  must  be  emphasized 
that  a  design  of  a  large,  complex  system  is  characterized  by  many 


design  changes  through  iterations.  As  a  result,  decision  support 
for  design  decomposition  should  accommodate  these  design 
alterations  brought  about  by  new  information  passed  between 
levels  of  the  network. 

1 . 1  The  Rationale  Behind  Decision  Support 

Decision  support  for  the  design  process  has  several 
appealing  functions  besides  potentially  improving  the  ultimate 
system  design.  A  decision  support  system  (  DSS  «•'  'uld  provide 
greater  managerial  control  of  the  design  proces  example,  a 
DSS  should  be  a  mechanism  for  coordinating  sut  /s  ;m  designs. 
This  would  allow  the  designer  to  recogni2e  when  a  particular 
subsystem  design  is  insufficient  to  meet  system  needs,  and  what 
the  best  alternative  designs  are.  This  feature  would  generate 
critical  information  early  in  the  planning  phase,  and  therefore 
direct  management’s  attention  to  the  important  design  issues. 
In  addition,  using  a  support  system  would  allow  the  designer  to 
determine  the  system  degradation  due  to  not  meeting  a  sped  lie 
system  specification. 

The  purpose  of  this  paper  is  to  describe  a  goal 
decomposition  approach  that  can  be  useful  in  a  DSS  for  the 
planning  phase  of  engineering  design.  In  the  following  section 
we  describe  some  of  the  approaches  that  currently  exist  for 
decision  support  of  the  design  process.  In  section  3,  a  goo i 
decomposition  method  is  examined  and  related  to  the  network  in 
Figure  1.  In  addition,  some  refinements  to  the  goal 
decomposition  algorithm  are  suggested.  Section  4  contains  two 
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example  applications.  Section  5  contains  the  concluding  remarks. 

2.  decision  Support  for  the  design  thocess 

The  engineering  design  problem  is  typically  expressed  as: 

(Pi)  min  f ( x , p ) 

subject  to:  g^x.p)  <0,  i  =  1 ,  .  .  .  ,  I  , 

h  •  ( x , b )  =  0,  j  =  1,...,J. 

J 

In  this  model,  the  vector  x  is  a  vector  of  design  variables;  that 
is,  these  are  the  variables  that  represent  the  freedoms  of  choice 
that  the  design  engineer  has  available.  The  vector  p  is  a  vector 
of  model  parameters. 

Johnson  and  Benson  [8],  (9J,  have  decomposed  the  design 

problem  by  treating  a  subset  of  variables  in  the  design  problem 
as  parameters  (g)  and  optimizing  over  the  remaining  variables. 
The  original  subset  of  "parameters"  are  then  released  and 
optimized. 

SobieSzczanski'Sobieski,  e  t  .  a  1 .  ,  [13J,  [14],  have  used 

linear  fjevompvs  i  t  j  on  in  an  attempt  to  simplify  the  large 
nonlinear  models  that  frequently  represent  a  design  problem. 
This  decomposition  is  accompl  islied  by  computing  the  sensitivity 
of  the*  global  solution  to  changes  in  sub  problem  solutions  (lf)]t 
and  representing  these  changes  by  first  order  approximations. 
The  ultimate  advantages  of  this  approach  result  from  the  explicit 
inclusion  of  oubproblem  constraints  in  the  global  optimization 
model  and  a  reduction  in  the  number  of  variables  in  the  global 
model.  in  recent  work  [ 6 ] ,  H3ftka  has  attempted  to  improve  upon 
the  linear  decomposition  method  by  removing  the  difficulties 
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encountered  due  to  discontinuities  in  the  derivatives  that 


migrate  upwards  through  the  hierarchy. 

Several  authors  have  suggest,  the  use  of  knowledge-based 
systems  to  aid  in  the  design  construction.  These  approaches  can 
be  dichotomized  by  the  type  of  knowledge  required  in  the  system 
and  the  use  of  that  knowledge.  In  one  case,  the  knowledge  is 
used  for  optimization  [1],  [2J.  Systems  of  this  type  can  be 
thought  of  as  addressing  the  more  general  issue  of  how  knowledge 
can  be  used  to  aid  in  the  search  of  optimal  solutions  for 
mathematical  models  that  are  characterized  by  non  1 i near  i  t 1 es . 
Surprisingly,  the  use  of  knowledge  specific  to  the  design  domain 
has  not  been  used  to  full  advantage  in  these  systems. 

The  other  type  of  knowledge-based  system  contains  knowledge 
specific  to  design,  but  does  not  emphasize  optimization  [ 3  j  ,  [ 4  j  . 
While  these  systems  can  produce  quality  designs,  the  designs  are 
based  on  heuristics  and  are  therefore  best  suited  for  producing 
designs  that  meet  stated  specifications. 

3.  GOAL  DECOMPOSITION 

In  this  paper  we  will  examine  a  method  of  decomposition  that 
makes  use  of  goal  programming.  One  of  the  primary  motivating 
factors  for  this  approach  is  that  it  is  felt  that  the  employment, 
of  goals  instead  of  "pure  constraints"  conforms  closely  to  the 
type  of  decision  making  currently  used  in  the  initial  design 
process.  Moreover,  ti.s  decision  support  architecture  which  is 
examined  may  take  advantage  of  a  variety  of  other  model  tng 
techniques  currently  available;  i.e.,  the  Sobieszczanski-Sobieski 
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approach,  knowledge-based  models,  etc. 

In  what,  follows,  we  will  use  the  term  objective  to  mean  a 
state  tli at  a  decision  maker  lias  identified  as  desirable  to 
attain.  The  term  attribute  is  used  to  describe  a  measure  which 
can  be  used  to  determine  the  degree  to  which  the  desirable  state 
has  been  achieved.  We  will  assume  throughout  this  paper  that 
with  sufficient  perseverance,  a  set  of  attributes  may  be 
determined  that  will  measure  the  objective  to  the  satisfaction  of 
the  decision  maker. 

3 . I  Dsfi^ioD 

One  of  the  fundamental  characteristics  that  differentiates 
the  goal  decomposition  approach  from  the  methods  of 
decomposition  that  have  been  found  in  the  median i ca I / ae r os  pace 
engineering  literature  is  the  treatment  of  design  variables. 
Consider  once  again  the  design  hierarchy  depicted  in  Figure  1. 
Imagine  that  an  optimization  has  been  performed  at  level  1. 
Information  has  been  passed  down  from  level  1  to  level  2  via  the 
vector  u.  This  vector  contains  the  attributes  of  the 
subcomponent  being  designed  at  level  2  that  are  desired  from  the 
perspective  of  the  component  being  designed  at  level  1.  As  a 
result,  the  vector  u  describes  what  it  is  that  we  want  buiLt  (the 
subcomponent  at  level  2). 

The  optimization  at  level  2  in  I ^  11  describe  how  we 
will  build  that  subcomponent  by  fur.,  c  refinement  of  the 
subcomponent,  into  its  major  subcomponents.  The  design  variables 
at  level  2  are  a  vector  w  that  contain  the  desired  attributes  of 
the  subcomponent  at  the  3rd  level  of  the  hierarchy.  These  values 
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are  determined  from  an  optimization  that  will  choose  w  so  that 
the  attributes  of  Level  1  (u)  are  met  as  closely  as  possible 
while  not  violating  any  constraints  that  are  placed  on  w. 

The  goal  programming  formulation  at  level  i,  i  =  1 . I  - 1 , 

may  be  described  as 

( F  2 )  min  Z  +  b~d"  (1) 

such  that: 

fk(w)  -  d^  +  d^  =  uk>  k  -  1  ,  •  •  .  ,  K ,  (11) 

w  e  0  .  (  3  ) 


As  described  above,  the  value  of  u^.  in  (2)  was  determined  from  a 

goal  optimization  at  level  ( i  —  1 ) ,  i  =  2 . 1-1.  In  the  case 

that  the  value  of  was  generated  at  the  Oth  level  (the  top 

level),  some  other  rau 1 t i ob j ec t i ve  technique  might  have  been  used 
to  determine  the  desired  values  of  the  attributes  of  the  top 
level . 


The 

values  b . 

k 

and  b^  arc  weights 

assigned 

t  o 

t.  he 

dev i  a  t ions 

dk  and  d k ' 

respectively,  which 

are 

the 

am oun t 

by 

which  we 

have  exceeded  or  fallen  short  of 

the 

kth 

g  o  a  l . 

The 

function  f^(w)  is  the  link  between  the  subcomponent  attributes 
and  the  component  attribute  . 

The  actual  value  of  the  kth  attribute  tha;.  cun  be  achieved 
in  (P2)  is  the  value  of  f^(w)  in  equation  (2)  (recall  that  the 
actual  value  that  we  would  like  to  achieve  is  u  ^  *  )  .  The  value  of 
the  kth  attribute  achieved,  f^lw*),  where  w*  is  the  optimal  value 
of  w  in  (P2),  may  not  equal  u^*  for  two  reasons:  (1)  The  goal  set 
designated  by  equations  (2)  may  be  overly  demanding  and  therefore 
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all  of  the  goal  may  not  be  achievable;  (2)  The  space  0  in 
equation  (3)  may  constrain  the  problem  from  achieving  all  of  the 
goals. 

The  goal  decomposition  algorithm  can  be  stated  as  follows: 

Al.  Optimize  the  attributes  of  level  i,  i  =  1,...,I-1, 
by  mathematically  relating  those  attributes  to 
attributes  of  the  subcomponen ts  (level  ( i  —  1 ) )  . 

A2.  Using  the  values  of  the  level  i  attributes  obtained 
in  the  optimization  at  level  (i-1)  as  goals, 
optimize  the  attributes  of  level  (i+1), 
i  =  1 . 1-1. 

A3.  If  necessary,  perform  a  reoptimization  at  level 
i  by  passing  constraints  on  the  attributes 
from  level  (i+1)  to  level  i. 

Several  iterations  between  levels  (i+1)  and  i  may  be  necessary  to 
obtain  a  satisfactory  design.  The  constraints  that  are  added  to 
the  level  i  model  in  step  A3  are  needed  to  cause  a  reallocation 
of  resources.  An  example  of  the  addition  of  these  constraints  is 
shown  in  the  first  application  of  Section  4. 

Since  all  of  the  goals  in  (F2)  may  not  be  met,  the  weighting 
and  in  the  objective  function  will  tend  to  dictate  which 
goals  are  satisfied.  Next  we  will  show  how  these  weighting  may 
be  selected  so  that  the  pursuit  of  the  goals  at  level  i  is 
consistent  with  the  desired  design  attributes  at  level  <  i-1). 

Theorem  1:  Let  (P2)  describe  tie  optimization  t  level  i  and  let. 

the  following  model  describe  the  descendant  optimization  at  levil 
(  i  +  l)  : 


R 


( P  3 ) 

such  that ; 


_  +  +  -  - 

min  £  c.e.  +  c.e. 

J  J  J  J 


g  ( x  )  -  e  .  +  e  . 
J  J  J 


=  w  .  ,  j  =  l . J, 

J 

x  e  X  . 


Then,  to  first  order  approximation,  the  values  of  c.  and  c 


should  be  chosen  according  to  the  following  algorithm: 

(Bl)  For  each  f  k  (  w  *  )  in  equation  (2),  use  th< 


compute 

new 

cons  t  an  t  s 

+ 

ak 

and 

ak 

+ 

+  -  - 

(Bl  .  1) 

If 

fk(w*)  = 

V 

let 

ak 

+ 

=  bk,  ak  =  bk ; 

+  -  + 

(  B  1  .  2  ) 

If 

fk(B*>  > 

uk’ 

let 

ak 

+ 

=  bk’  ak  =  'bk; 

(  B  1  .  3  ) 

If 

fk(w*)  < 

V 

let 

ak 

=  -bk,  ak  =  bk. 

Form 

the 

values  c+ 

J 

and 

c  . 

J 

i  n 

equation  (4)  by 

c+  =  £  3f,/3w.  (w*)  , 

j  k  k  j 


c.  =  E  a,  5f./aw.  (w*). 


Proof : 


foil ow  i  ng 


Consider  the  first  order  Taylor’s  series  expansion  of  f^(w) 


around  w*  : 


fc(w)  “  fu(Sf*) 


E  4f./3w  (w*)  (w  -  w  .  ) 

**  J  J  J 


The  term  on  tlie  left  hand  side  of  (9)  is  the  amount  that  the 

optimal  value  f^(w*)  has  been  missed  by  the  selection  of  w.  The 

right  hand  side  of  (9)  contains  the  term  ( w  .  -  w.  ) ,  which  is 

■J  J 

the  amount  by  which  the  Level  (i+l)  optimization  was  unable  to 


meet  the  goal  w^.  that  was  passed  down  from  the  level  i 
optimization.  Thus  equation  (9)  depicts  the  influence  that  the 
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failure  to  meet  goal  w*  at  level  (i+1)  has  on  the  optimization 

J 

♦ 

at  level  i.  If,  for  example,  f^(w*)  >  u^  ,  the  optimal  value  of 
w*  caused  the  attribute  represented  by  f^(w*)  at  level  i  to 
be  exceeded,  and  d£  >  0,  d^  =  0  in  equation  (2).  If  the  value 
of  w  realized  from 

w  .  =  g  .  (x) ,  j  =  1 . J, 

J  v/ 

does  not  equal  w  (the  goal)  but  instead  exceeds  w.  ,  then  the 

J  d 

value  of  the  objective  function  at  level  i  will  change  by 


bk  ffk(y)  “  fk(y*)j* 

(here  we  assume  that  w  is  near  w*  so  that  f^(w)  >  u^*)* 

But  this  is  just  the  right  hand  side  of  (9)  multiplied  by  the 
appropriate  weight  from  the  level  i  optimization.  Multiplying 
the  right  hand  side  of  (9)  by  b*  yields  the  result  for  the 
special  case  that  f^(w*)  >  u^*  (number  2,  part  A  of  the 

algorithm).  The  other  parts  of  the  algorithm  can  be  determined 
from  similar  arguments. 


3.2  Refinements  to  the  Goal  Decomposition  Algorithm 

(1)  If  there  are  no  constraints  in  equation  (3)  that 

connect  w.’s  from  different  subproblems,  then  a  feasible  design 
J 

at  level  i  may  be  assured  easily  at  level  (  i  *■  1  )  by  adding  the 
constraint 

g  j  (  x )  e  0 
to  problem  ( P3 ) . 

(2)  If  the  optimization  at  level  (i+1)  does  not  meet  the 
goals  w*,  then  a  reoptimization  at  level  l  may  generate  a  better 
(or  feasible)  solution  when  taking  into  account  the  best.  valuer. 
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of  the  attributes  that  could  be  achieved  at  level  (i+1). 


In  this 


following  cuts 

may  be 

added 

to 

equation 

(A) 

If 

8j(x.) 

=  w  .*  , 

J 

then 

add 

no  cuts ; 

(R  ) 

I  f 

8j(x.) 

>  w  .  *  , 

J 

t.ien 

add 

the  cut 

w  .  > 

J 

; 

(C  ) 

If 

gj(x«) 

<  w  .  *  , 

J 

then 

add 

the  cu  L 

w  .  < 

J 

8j<x*> 

• 

(3)  The  problem  that  a  multipurpose  component  (a  component 
with  multiple  parent  nodes  in  the  decomposition  network)  causes 
can  be  coordinated  in  this  process  by  adding  separate  goals  from 
each  parent  optimization.  For  example,  if  a  multipurpose 

component  is  passed  down  the  goal  w  ^  *  =  1  from  one  parent 

optimization  and  w^*  =  2  from  the  other  parent  optimization, 
then  the  following  two  goals  should  be  added  to  the  subproblem: 

g^x)  +  c4  ~dl  =  1  • 
g^Cx)  +  d2  ~  d2  =  2> 

where  w j  =  gj(x).  The  coefficients  in  the  objective  function 
will  help  determine  which  of  these  contradicting  goals  should 
dominate. 


4.  applications 

4.1  A  Conceptual  Application 

The  first  system  design  that  will  be  optimized  using  goal 
decomposition  is  depicted  in  Figure  2a.  We  begin  with  this 
conceptual  system  to  illustrate  the  techniques  that  may  be 
employed  to  model  system  cost  and  reliability.  Similar  models 
may  bo  built  to  include  attributes  measuring  system  performance, 


supportability,  and  producibiLity. 

The  network  decomposition  of  this  system  is  shown  in  Figure 

2b.  At  level  1,  the  system  is  viewed  as  consisting  of  two 

components  in  series,  as  indicated  by  the  dotted  line  in  Figure 

2a.  We  will  assume  that  Component  1  may  be  viewed  as  consisting 

of  two  identical  components  in  parallel,  while  component  2 

consists  of  two  components  in  series. 

We  shall  also  assume  that  an  optimization  has  been  performed 

at  level  0  that  specified  the  target  values  of  the  reliability 

$  $ 

and  cost  of  the  entire  system  to  be  R  =  .995  and  C  =  $2000. 

Let: 

=  Reliability  of  component  i,  i  =  1,2, 

=  Probability  that  the  component  operates  without 
failure  over  a  specified  time  interval, 

=  Cost  of  component  i,  i  =  1,2, 

R.  .  =  Reliability  of  subcomponent  (i,j),  i  =  1,2, 

1J  j  =  1,2, 

C  .  .  =  Cost  of  subcomponent  (i,j),  i  =  1,2,  j  =  1,2. 

■i  j 

It  is  also  assumed  that  Rj  and  roust  be  at  least  .95  and  .1)2 
respectively,  and  that  the  total  budget  for  the  system  is  $2600. 

In  order  to  build  a  mathematical  model  of  the  system,  it  is 
first  necessary  to  determine  how  reliability  influences  cost. 
Suppose  that  by  experience  it  is  known  that  the  reliability  of 
the  subcomponents  is  related  to  the  cost  of  the  subcomponents  in 
the  following  way: 

Cn  -  -200  In  {  l  -  Rh(2  -  Rn)}  ,  (10) 

C  jo  =  C  j  j  ,  ill) 

C21  =  1 . 2/ ( 1  -  H21) ,  (12) 
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(13) 


C22  ~  1/^1  R22)‘ 

Notice  that  in  each  of  these  cases  the  cost  of  the  subcomponent 

increases  rapidly  as  the  reliability  approaches  1. 

The  development  of  the  cost  functions  and  C9  for  an 

optimization  at  level  1  in  Figure  2b  presents  two  significantly 

different  problems.  Component  1  consists  of  two  identical 

components  in  parallel;  and,  therefore,  the  cost  function  may  be 

written  simply  as 


since 


'  C11  f  C12 


=  -400  1 n { 1  -  R  }  , 


(14) 


R1  "  Rll(2  Rll) ’ 

On  the  other  hand, 

C2  =  C21  *  C22 

cannot  be  written  as  a  simple  function  of 


( in: 


(lb) 


R2  '  R2 1R22  ‘  (  l7) 
Consequently,  the  developer  of  a  model  has  two  alternative 

approaches;  (l)  Include  in  the  level  1  optimization  model 

mathematical  expressions  that  involve  and  R ^ ,  which  are 

level  2  attributes;  or  (2)  Estimate  the  cost  function  Cr>  as  a 

function  of  -  The  former  alternative  regresses  toward  the 

problem  that  prompted  decomposition  in  the  first  place:  the 

model  begins  to  grow  larger  because  we  have  decided  to  include 

too  much  detail  in  the  model  at  the  higher  levels  in  the  network. 

The  second  alternative  risks  the  loss  of  some  information  in  the 

estimation  process  but.  maintains  the  spirit,  of  decomposition. 

For  example,  using  equations  (12)  and  vl3),  the  following 
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table  of  costs  for  alternative  designs  may  be  easily  constructed: 


R2i  (1  = 

1,2)  C21 

C22 

.  9 

$  12 

$  10 

.  925 

16 

13.3 

.  95 

24 

20 

.  975 

48 

40 

.  99 

120 

100 

.  995 

240 

200 

.  999 

1200 

1000 

Table  I'-Cg^  anc*  ^°2 


Using  equation  (17)  and  different  combinations  of  R , 


ami 


R 22 >  values  of  may  be  computed  from  Table  1.  Next,  a  least 
square  estimation  technique  may  be  used  to  develop  an  estimate;  of 
C9  as  a  function  of  R2 .  Such  a  model  was  developed  for  t h i s 
example  from  15  combinations  of  and  R ^  which  yielded  the 


cost  function 


=  154.2  +  4.8(1  -  R_) 


-1 


(  18 


The  optimization  model  at  level  1  may  now  be  written  as: 
(Ml )  min  2 

such  that 


5000di  +  d2 


K1R2  “ 
C 1  +  C2 
C1  +  C2 
*\d\ 
d2d2 


dl  +  dl 


*2 


d  „  +  d ,, 


=  .  995 

=  2000 
£  2600 
=  0 


.95  <  R  £  1 
.92  £  R.,  £  1 


)  4 


0 


as  well  as  equations  (14)  and  (18)  are  satisfied  and  all  of  the 
variables  in  the  model  are  non-negative. 

The  coefficients  in  the  objective  function  of  (Ml)  reflect 
the  trade-ofis  that  were  judged  reasonable  by  the  decision  makers 
at  level  0  and  have  been  assumed  for  this  example.  The  optimal 
solution  to  (Ml)  generated  by  GINO  [10]  is: 

z  =  158.187  d"  =  °  C  2  =  455.48 
d~  =  .082  d^  =  0  =  .97896 
dj  =  0  C  =  1544.52  J?2  =  .98407 

4.1.1  Level  2  Optimization  Models 

To  generate  a  model  that  optimizes  the  design  of  the 
subcomponents,  we  need  to  use  the  values  of  and  R^,  i  ],2, 
that  were  computed  in  the  optimal  solution  to  (Ml).  In  addition, 
the  results  of  Theorem  1  are  needed  to  search  for  solutions  thal 
are  consistent  with  the  level  1  optimization.  The  models  (Mil) 
and  (MI2)  for  components  1  and  2  respectively  and  their 
solutions  may  be  found  in  appendix  A. 

(Mil)  and  (M12)  contain  constraints  that  are  local  to  the 
subcomponent  design  as  depicted  by  equation  (6).  For  example,  in 
(Mil),  l?n  <  .8.  The  solutions  for  C  ^  and  R  ^  in  (Ml)  are  used  as 
the  goals  in  (Mil)  and  (M12).  Also  notice  that  the  constraints 
0.95  <  R^  <  1  have  been  included  in  (Mil)  as  suggested  in  (l)  of 
section  3.2. 

Table  2  compares  the  values  of  FL  and  C^,  i  -  1,2,  t. hut  were 
generated  in  the  solutions  of  (Ml)  with  (Mil)  and  ( M 1 2 ) .  It  is 
clear  from  this  table  that  the  largest  value  of  that  can 
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possibly  be  generated  is  .96,  and  that  $1267.55  is  the  largest 
budget  needed;  otherwise,  additional  dollars  would  have  caused 
the  solution  of  R  ^  in  (Mil)  to  be  larger. 


Variable 

De  r i ved  f  rom 

Level  1 

Achieved  at 
Level  2 

Rl 

. 97896 

.  96 

C1 

1544 .52 

1287 . 55 

R2 

. 98407 

. 97752 

C  2 

455.48 

455.48 

Table  2--Attribute  Goals  and  Values  Realized 
in  (Mil)  and  (M12) 


If  the  constraints  R^  <  .96  and  <  1287.55  are  added  to 
(Ml),  the  new  solution  to  (Ml)  becomes 


z  = 

2 16. 272 

d2 

=  0 

C2  =  712.45 

dI 

=  .  043 

d.t 

=  0 

R2  =  .96 

dl 

=  0 

C1 

=  1287.55 

R2  =  .9914 

This  solution  is,  in  fact,  the  globally  optimal  design;  i.e.,  the 
design  that  would  be  generated  if  we  had  built  one  model  that  had 
contained  all  of  the  constraints  on  the  component  and 
subcomponent  designs. 

In  summary,  in  order  to  generate  the  globally  optimal 
design,  it  was  necessary  to  pass  information  oownward  through  the 
decomposition  network  that  described  the  desired  attributes  of 
the  system.  Information  in  the  form  of  constraints  was  passed 
upward  through  the  network  that  indicated  the  best  achievable 
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values  of  the  attributes.  A  reoptimization  was  performed  that 
reallocated  the  resources  (in  this  case  budget  dollars)  to 
improve  the  design.  Further  optimization  iterations,  in  general, 
would  continue  to  improve  the  design. 

4.2  A  Second  Application 

Consider  the  problem  of  designing  a  water  pumping  system  and 
a  structural  support  for  the  pumping  system.  As  shown  in  Figure 
3a,  it  is  assumed  that  the  water  pumping  system  consists  of  two 
parallel  pumps  that  draw  water  from  a  lower  reservoir  and  pump 
the  water  into  a  reservoir  that  is  150  feet  above  the  lower 
reservoir  [16].  The  pumping  system  is  to  be  positioned  on  a  30 
foot  reinforced  concrete  beam  of  depth  30  inches  as  discussed  in 
[9].  The  decomposition  network  for  the  pump/beara  system  is  shown 
in  Figure  3b. 

Assume  that  a  budget  of  $950  may  be  spent  on  the  project  but 
that  the  desired  budget  is  $900  (as  determined  in  an  optimization 
at  level  0).  A  goal  of  200  feet  for  the  head  loss  has  been 
established  to  take  into  account  the  loss  of  energy  due  to 
friction  in  the  pipes.  This  head  loss  must  be,  at  a  minimum,  160 
feet.  In  addition,  a  goal  of  6480  in  k  of  flexural  strength  is 
desired  in  the  beam.  Define: 

Cj  -  Cost  of  the  pump  subsystem, 

C0  -  Cost  of  the  beam  subsystem, 

H  "  Head  loss, 

S  -  Flexural  strength  of  the  beam. 

Prior  experience  has  shown  that  C^  and  increase  with  H 
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and  S  in  the  following  way: 


C1  = 

.01125H2  +  . B75H , 

(19) 

Co  = 

. 05787S . 

(30) 

Assuming  once  again  that  the  optimization  at  level  0  has  supplied 
the  coefficients  in  the  objective  function  that  reflect  the 
trade-offs  that  would  be  considered  appropriate  at  the  top 
level,  the  purap/beam  optimization  model  can  be  written  as: 


(M2)  min  z  =  d*  +  5d2  +  . 0926d~ 

such  that 


C1  +  C2 

’  dl 

+  dl 

- 

900 

H 

+ 

-  do 

fa 

+  d2 

rr 

200 

S 

“  d3 

+  d3 

= 

6480 

C1  +  C2 

< 

950 

H 

> 

160 

d  +  d~ 

l  i 

= 

0,  i  =  1,2,3, 

equat  ions 

(19)  and  (20)  hold 

and  all  variables 

non negative.  The 

optimal  so 

lution  to  (M2)  is: 

z  = 

96.875 

■ 

0 

C,  ■  375 

di  = 

0 

d  3  = 

0 

H  =183.34 

<  - 

13.59 

d;  ■ 

0 

S  =  6480 

d„  =  16.66  C.  =  538.59 

fa  1 


in  the  model 


are 


4.2.1  Level  2  Models 


Next  consider  the  optimization  of  the  pump  subject  to  local 
constraints.  In  this  example,  assume  that  the  headless  may  be 
expressed  in  terms  of  the  flow  rates  (cubic  feet/se condj,  IV, 
through  pump  i ,  i  -  1,2: 
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H 

= 

A 

+  30  Fj 

-  6F2 

(21  ) 

H 

B 

+  20F2 

~  12F2 

(22) 

In  these 

expressions,  A 

and  B  are  parameters 

which 

specify 

attributes 

of 

the  pumps. 

The 

selection  of 

the 

pumps  is 

res  t  r icted 

t  o 

50  < 

A 

<  130 

(23) 

n 

4+ 

00  < 

B 

<  220. 

(24) 

The  amount 

of 

f  l 

ow  in  each 

pump  is 

restricted  to 

be  at 

least.  2 

cubic  feet/second. 


Let  Cjj  be  the  cost  of  pump  j,  j  =  1,2.  The  cost  of  the 

pumps  is  related  to  the  flow  rate  via 


C11  - 

100  r 

(25) 

C  1 2  s 

120  F2. 

(26; 

The 

total  budget 

for  the  pump 

subsystem  will 

be  at  least  1525. 

The 

optimization 

model  (M21) 

for  the  pump 

subsystem  and  the 

optimal  solution  is  given  in  appendix  B. 

The  beam  subsystem  has  design  variables  W  and  V,  which 
represent  the  width  of  the  beam  and  the  cross  sectional  area  cl 
the  reinforcing  steel.  For  this  example,  it  is  assumed  that  the 
reinforcing  steel  has  standard  sizes  that  come  in  .'15  inch 
increments.  The  flexural  strength,  S,  of  the  beam  is  related  to 
W  and  V  by 

V  -  .2458  V2/ W  >  S/i080.  (27;. 

The  unit  cost  of  the  concrete  is  $. 02/s qua  re  inch /lineal 
foot,  and  the  unit  cost  of  steel  is  $  1  .  00/square  inch  /I  in  cal  foot. 
The  length  of  the  reinforcing  steel  is  slightly  less  than  the 
length  of  the  beam  (29.4  feet),  so  that  the  cost  of  Lhe  steel  is 
given  by  29.4V.  The  cost  of  the  concrete  is  given  by 
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,02  x  (cross-sectional  area  of  the  beam)  x  30 
=  .02  x  (W  x  30)  x  30 

=  18  W. 

The  width  of  the  beam  must  be  at  least  8  inches  to  provide 
sufficient  support  to  the  pump  subsystem.  Once  again  using 
Theorem  l,  the  optimization  model  (M22)  for  the  beam  subsystem 
may  be  formulated  and  solved  as  shown  in  appendix  C. 

The  value  of  V  in  the  solution  of  (M22)  is  not  a  integral 
multiple  of  .25;  as  a  result,  this  solution  is  not  truly  feasible 
given  the  problem  constraints.  In  order  to  achieve  the  optimal 
feasible  solution,  two  new  models  (M23)  and  (M24)  may  be  formed 
from  (M2 2)  bv  adding  the  constraints  V  <  7.75  in  the  first  case 
and  V  >  8  in  the  second  case.  Model  (M24)  has  no  feasible 

solution,  while  the  optimal  solution  to  (M23)  is  given  by: 


7.  =12.69 

-2  -  3.15 

S 

=  6376.32 

+ 

eL  =  103.05 

e2  =  0 

W 

=  8 

“f 

el  =  ° 

C2  =  371 . 85 

V 

=  7.75 

This  solution  is  therefore  the  optimal  attributes  for  the  design 
of  the  beam  subsystem.  The  solution  to  (M21)  and  (M23)  may  be 
used  to  generate  constraints  in  order  to  further  constrain  the 
solution  to  (M2).  This  would  begin  a  new  iteration  of  the  goal 
decomposition  algorithm.  However,  a  designer  may  decide  that  t he 
design  achieved  by  combining  the  solution  to  i M2 1  )  and  (M23)  to 
he  satisfactory. 


5.  CONCLUDING  RKMARKS 

This  paper  has  attempted  to  develop  a  broader  perspective 
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of  the  engineering  design  problem  by  distinguishing  the  type  of 
optimization  problems  that  occur  at  the  top  and  lower  levels  of 
the  design  hierarchy.  One  key  feature  of  middle  and  lower  level 
optimizations  is  that  subsystem  optimization  does  not  correspond 
directly  to  performance,  cost,  etc.  optimization.  A  goal 
decomposition  method  was  described  and  illustrated  for 
subcomponent  optimization.  This  decomposition  method  had  a 
natural  objective  function  that  tended  to  cause  the  solution  of 
the  subproblem  to  drive  the  solution  of  the  optimization  one 
level  higher  to  its  best  possible  feasible  solution. 

The  architecture  of  the  decision  support  system  described 
using  goal  decomposition  readily  enhanced  the  iterative  nature  of 
the  design  process.  In  the  first  example  shown  in  this  paper,  a 
single  reoptimization  found  the  best  global  solution  to  the 
design  problem  once  additional  information  was  supplied  to  the 
level  1  model  by  subsystem  optimizations. 
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APPENDIX  A 


Optimal  solution: 
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R22 

> 

99 

R2  2 

>  1 

+  - 
e  .  e  . 

1  3 

=  0 

i 

=  1.2 

as  well 

as  equations 

(  12)  , 

(13),  and 

(  16) 

and  all  variable 

non  n  etf a  t  ive. 

Optimal 

solution: 

2 

=  32.04 

+ 

e2 

=  0 

R22 

=  .99747 
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=  0 
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APPENDIX  B 


( M2  I )  rain 
such  that 


z  = 


-5e^  +  5e^  + 


H  "  el  +  el 

. 98407 

C1  "  e2  +  6  2 
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538.59 

C1 
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Cli  +  C12 

C1 
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525 
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F2 
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+  - 
e  .  e  . 

i  l 

= 

0  i  =  1 , 2 

as  we  1  1 

as  equations  (  2  l  )  -  (  2  6  ) 

and  all 

variables  are 

nonnegat ive 

Opt ima  1 

solution: 

z 

=  73.25  C{  = 

525 

F  =4.71 

+ 

e  1 

s  °  C11  = 
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F1  =  2 

el 

=  17.37  Cl2  = 
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F2  =  2.71 
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e2 

=  0  A 

130 

H  =165. 
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-  13.59  B 
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APPENDIX  C 


(M22) 

min  z  =  .0926e^  - 

s  uch 

that 

S  ~  e  j  +  e  j 

=  6480 

C2  "  e2  +  e2 

=  375 

C2 

=  29.4V  +  18W 

W 

|\z* 

00 

+  - 
e  .  e 
i  i 

-  0  i  =  1,2 

as  well  as  equation  (27)  and  all  variables  are  nonnegative. 


Optimal  solution: 
z  =  12.69 

el  =  0 
ex  =  103.05 

e0  =  0 

U 

e  t)  -  3.15 


C0  =  375 
S  =  6437.17 
W  =  8 
V  -  7.857 
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Figure  1 — Decomposition  Network 
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Figure  2a — A  Conceptual  System 
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Figure  3b — Network  Decomposition  of  the 
Water  Pumping  Systsn 


